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Abstract. Starting from Holstein's work on small polaron hopping, the evolution equations for localized 
and extended states in the presence of atomic vibrations are systematically derived for an amorphous 
semiconductor. The transition probabilities are obtained for transitions between all combinations of local- 
ized and extended states. For any transition process involving a localized state, the activation energy is 
not simply the energy difference between the final and initial states; the reorganization energy of atomic 
configuration is also included as an important part of the activation energy (Marcus form). The activa- 
tion energy for the transitions between localized states decreases with rising temperature and leads to the 
Meyer-Neldel rule. The predicted Meyer-Neldel temperatures are consistent with observations in several 
materials. The computed field-dependence of conductivity agrees with experimental data. The present 
work suggests that the upper temperature limit of variable range hopping is proportional to the frequency 
of first peak of phonon spectrum. We have also improved the description of the photocurrent decay at low 
temperatures. Analysis of the transition probability from an extended state to a localized state suggests 
that there exists a short-lifetime belt of extended states inside the conduction band or valence band. 

Key words, polaron - short-lifetime belt - electron transfer - reorganization energy - conductivity 

PACS. 72.15.Rn Localization effects, 72.20.Ee hopping transport, 73.20.Jc Delocalization processes, 
72.80.Ng Disordered solids. 



1 Introduction 

■ In the last 50 years, transport properties in amorphous 
semiconductors have been intensely studied [IJ[2j[3j[4j[5l[6] . 
Miller- Abrahams (MA) theory[7] and variable range hop- 
ping (VRH) [S] are frequently used to fit dc conductiv- 
ity data. For localized tail states close to a mobility edge, 
'phonon induced delocalization' [9, plays an important role 
in transport. In addition, exciton hopping among local- 
ized states is suggested as the mechanism of photolumi- 
nescence in a quantum well [10 . The Meyer-Neldel rule has 
been deduced from the shift of Fermi level |llj. from multi- 
excitation entropy [T2I[T5] and from other perspectives [T31 

; using. 

In the MA theory of dc conductivity [7], the polar- 
ization of the network by impurity atoms and by carri- 
ers in localized tail states was neglected. The transitions 
between localized states were induced by singlc-phonon 
absorption or emission. Subsequent research on transient 
photocurrent decav [T71ll"Blll9ll20] adopted a parameterized 
MA transition probability and therefore inherited the single- 
phonon features of MA theory. However in other elec- 
tronic hopping processes, the polarization of the environ- 
ment by moving electrons plays an important role. Elec- 
tron transfer in polar solvents, electron transfer inside 



large molecules [H] and polaron diffusion in a molecular 
crystal 22, 23] are relevant examples. The transition prob- 
ability Wll between two sites in a thermally activated 
process is given by the Marcus formula [2Tj: 

W LL = VLL e- B S L "'* T , E^=^(l + ^f, (1) 

4 All 

where vll has the dimension of frequency that character- 
izes a specific hopping process. Here, E^ L is the temperature- 
dependent activation energy, All is the reorganization 
energy, AG\ L is the energy difference between the fi- 
nal state and the initial state[2T|. ((T|) has been estab- 
lished for both electron transfer [2"T] and small polaron 
hopping [24l. The mathematical form of Holstein's work 
for Id molecular crystals is quite flexible and can be used 
for 3d materials with slight modifications [2"51I26[I27II28II29| . 
Emin applied small polaron theory to transport proper- 
ties in amorphous semiconductors, and assumed that the 
static displacements of atoms induced by electron-phonon 
(e-ph) interaction caused carrier self-trapping on a sin- 
gle atom [3D]. The effect of static disorder was taken into 
account by replacing a fixed transfer integral with a dis- 
tribution of. He found that the static disorder reduces the 
strength of the electron- lattice coupling needed to stabilize 
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global small-polaron fomation[31 . Bottger and Bryksin 
summarized [32] a wide range literature on hopping con- 
duction in solids before 1984. 

By estimating the size of various contributions, we 
show in Sec. II that (1) static disorder is more important 
than the static displacement induced by e-ph interaction, 
so that the carriers are localized at the band tails by 
the static disorder pi3"3"j; (2) the static displacement in- 
duced by the e-ph interaction in a localized state is much 
larger than the vibrational amplitude of the atoms, so 
that any transition involving localized state(s) must be 
a multi-phonon process. In a semiconductor, the mid-gap 
states and localized band tail states are the low-lying ex- 
cited states (most important for transport) at moderate 
temperature[8 ( 33 . To simplify the problem, let us leave 
aside the mid-gap states (induced by impurity atoms, dan- 
gling bonds and other defects) and restrict attention to 
localized tail states (induced by topological disorder) and 
extended states only. Mid-gap states are well treated by 
other methods [5]. 

The first aim of this work is to extend Holstein's work[22 , 
23J to amorphous semiconductors, and to derive the equa- 
tions of time evolution for localized tail states and ex- 
tended states in the presence of atomic vibrations. If there 
are only two localized states in system, the evolution equa- 
tions reduce to the Marcus theory of electron transfer. If 
there is one localized state and one extended state, the 
evolution equations are simplified to Kramers' problem 
of particle escape-capture. The short-time solution of the 
evolution equations can be used to compute a spatially av- 
eraged current density [34 i.e. conductivity (we will report 
this in a forthcoming paper). 

The second aim of this paper is to estimate the tran- 
sition probabilities of four elementary processes: (i) tran- 
sition from a localized tail state to another localized tail 
state (LL); (ii) transition from a localized tail state to 
an extended state (LE); (iii) transition from an extended 
state to a localized tail state (EL) and (iv) transition from 
an extended state to another extended state (EE). The 
distribution functions of carriers in localized states and 
in extended states satisfy two coupled generalized Boltz- 
mann equations. The transition probabilities of LL, LE, 
EL and EE transitions obtained are necessary input for 
these two equations. Then transport properties could be 
computed from these Boltzmann equations. In this paper, 
we will not pursue this approach, and instead only esti- 
mate conductivity from the intuitive picture of hopping 
and mean free path. The present work on four transitions 
illuminates the physical processes in transport which are 
obscured in ab initio estimations of the conductivity [M]. 

The third aim is to (i) check if the new results describe 
experiments better than previous theories; (ii) establish 
new relations for existing data and (iii) predict new ob- 
servable phenomenon. In LL, LE and EL transitions, the 
frequency of first peak v of the phonon spectrum supplies 
energy separating different transport behavior. In the high 
temperature regime (T > hP/ks), static displacements in- 
duced by the e-ph interaction require configurational re- 
organization in a transition involving localized state (s). A 



Marcus type transition rate is found at 'very' high tem- 
perature (T > l.hhv /ks). The decreasing reorganization 
energy with rising temperature leads to the Meyer-Neldel 
rule: this is a special dynamic realization of the multi- 
excitation entropy model[T3]. The predicted Meyer-Neldel 
temperature is satisfied in various materials (cf. Table [1]). 
At low temperature (T < hP/lOks), the atomic static dis- 
placements induced by e-ph interaction reduce the transfer 
integral, and a few phonons may be involved in LL, LE 
and EL transitions. For LL transition: VRH may be more 
effective than the transition between neighboring local- 
ized states. The upper temperature limit of VRH is found 
to be proportional to P. This new relation is compared 
to existing experimental data in Figure [31 In the inter- 
mediate range(/i^/10fcs < T < hP/ks), the well-known 
non-Arrhenius and non-VRH behavior appears naturally 
in the present framework. Comparing with previous mod- 
els, the present work improves the description of the field- 
dependent conductivity, cf. Figure [1] and Figure [5J The 
predicted time decay of the photocurrent is compared with 
observations in a-Si:H and a-As2Se3 and is improved, com- 
pared to previous models at low temperature, cf. Figure 
U] For EL transition, we suggest that there exists a short- 
lifetime belt of extended states inside conduction band or 
valence band (cf. Figure [5]). 



2 Evolution of states driven by vibrations 

2.1 Single-electron approximation and evolution 
equations 

For definiteness, we consider electrons in the conduction 
band of an amorphous semiconductor. For carriers in mid- 
gap states and holes in the valence band, we need only 
modify the notation slightly. For the case of intrinsic and 
lightly doped n-type semiconductors, the number of elec- 
trons is much smaller than the number of localized states. 
The correlation between electrons in a hopping process 
and the screening caused by these electrons may be ne- 
glected. Essentially we have a single particle problem: one 
electron interacts with localized tail states and extended 
states in the conduction band. 

Consider then, one electron in an amorphous solid with 
Af atoms. For temperature T well below the melting point 
T m , the atoms execute small harmonic oscillations around 
their equilibrium positions {lZ n }: W n = lZ n + u n , where 
{W n } and {u 11 } are the instantaneous positions and vi- 
brational displacements. The wave function ip(r, {u n };t) 
is then a function of the vibrational displacements of all 
atoms and of the coordinates r of the electron. The Hamil- 
tonian H\ of the "one electron + many nuclei" system may 
be separated into: H\ = h e + h v , where: 

h e = —V 2 + J2u(r,K n ,u"), (2) 
2m 

n=l 

is the single-electron Hamiltonian including the vibrations. 
In h e put {u 11 = 0}, and one obtains the hamiltonian h a 



Ming-Liang Zhang and D. A. Drabold: Phonon driven transport in amorphous semiconductors 



3 



for an electron moving in a network with static disorder. 
Let 

h 2 1 

K = Y ~^" V j + 2 H k JkXjX k , (3) 
3 3 jk 

be the vibrational Hamiltonian, where (kjk) is the ma- 
trix of force constants. To simplify, we rename {u 11 , n = 
1, 2, • • • ,W} as {xj, j = 1, 2, • • • , 3W}- The evolution of 
tjj(r, {u n }; t) is given by: 

^ ^ (r 'f ; ^ Mr,K};t). (4) 

The Hilbert space of h e is spanned by the union of 
localized states {<^4i} and the extended states {£bi}- 
ip(r, {u 11 }; t) may be expanded as 

V>(r,aii, • • ■ ,x 3A r;t) = Y a M ( t>A 1 + ^&b 1 £b 1 , (5) 

where is the probability amplitude at moment t that 
the electron is in localized state Ai while the displace- 
ments of the nuclei are {xj, j — 1,2, •• • 3A/"}; bs 1 is the 
amplitude at moment t that the electron is in extended 
state B± while the displacements of the nuclei are {xj, j = 
1, 2, • • ■ 3Af}. To get the evolution equations for a Al and 
6_Bj, one substitutes (O into and separately applies 
/ d 3 r<^ 2 and J d 3 r£ B2 to both sides of the equation|22, 
[2"3"j . After some approximations (cf. Appendix A), the evo- 
lution equations are simplified to: 

d 

i^gl _ E M - h v}dA 2 = Y J MA 1 dAi + Y K A 2 B 1 b B 1 , 

Ai Si 

(6) 

and 

d 

i^gl ~ Eb 2 ~ h>v}bB 2 = Y, J / B 2 A 1 a A 1 + Y, ^B 2 B^Bii 

(7) 

where £?b 2 is the energy eigenvalue of /i a corresponding to 
extended state £b 2 , and 

E A 1 ({Xp 1 Ai }) = E° Ai - Yl d *±A ( 8 ) 
pa 1 eD Al 

is the energy of localized state (f>A 1 to the first order of e-ph 
interaction, where d PA = J dr\<j> Ai (r, {lZ n })\ 2 dU/dX PAi . 
E A and ^P A are the corresponding eigenvalue and eigen- 
function of h a and <j> Al is correction of <fP A to the first 
order of e-ph interaction. 
The transfer integral 

Ja 2Ai = j dM\ 2 Y U(r-K p ,u*>)<p Al , (9) 

induces transitions from <p Al to $a 2 - Similarly J AlA _ 2 in- 
duces transitions from A 2 to A\, and is due to the attrac- 
tion on the electron by the atoms in D Al . From definition 
([H]), no simple relation exists between J A2Al and J AlA2 . 



This is because: (i) the number of atoms in D Al may be 
different to that in D A2 ; (ii) even the numbers of atoms 
are the same, the atomic configurations can be different 
due to the topological disorder in amorphous materials. 
This is in contrast with the situation of small polarons in 
a crystal [25]: J AciAl — (Ja 1 A 2 )*> where two lattice sites 
are identical. The non-Hermiticity of Ja 2 Ai comes from 
the asymmetric potential energy partition (|59p for well lo- 
calized states. The transition probability from a localized 
state to another is different with its inverse process. 
The e-ph interaction 




is a linear function of atomic displacements Xj . It causes 
EL transitions from extended states to localized states. 
LE transition from a localized tail state in region D Al to 
an extended state (LE) is induced by the transfer integral: 

j 'b 2 a^ E [<Pr& a U(r-H p ,u?)<j>A 1 , (11) 

not by the e-ph interaction K' A B . Later we neglect the 
dependence of J'b 2 Ai 011 the displacements of atoms and 
only view J'b 2 a as function of £ Al only. The EE transition 
between two extended states £b x and £b 2 is caused by e-ph 
interaction: 




K B2Bl = (K BlB2 )* . (12) 

It is almost identical to scattering between two Bloch 
states by the e-ph interaction. In contrast with LL, LE 
and EL transitions, transition £g 1 — > ^b 2 and its inverse 
process £b 2 £bi are coupled by the same interaction, 
as illustrated in (TT2"1) . The transition probabilities of the 
two processes are equal. The numerical magnitude of the 
coupling parameters of the four transitions are estimated 
in Appendix B. 

2.2 Reformulation using Normal Coordinates 

As usual, it is convenient to convert {x k } to normal coordinates |37l 

Xk = Y ^kaOa, {A T kA)p a = SafjMaLol, (13) 
a 

a = 1,2,---,3A/", 

where A ka is the minor of the determinant | kit — u> 2 MiSik | = 
0, A T is the transpose matrix of the matrix (Ak a ). The 
two coupling constants in (ITU1) and (|12l) which involve e-ph 
interaction are expressed as: 

K 'a 2Bi =Y & <* K Zb 1 , (14) 
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and 

K B 2 Bx = ^ @ol K %B, » 
a 

Xg 2Bl = J2 A d° I d3 <k , (15) 

3 J 

where K'g B and K^ 2Bi ha,ve the dimension of force. ((SJ) 
and (O become: 

d 

(ih— - h A . 2 )a A2 (■■■O a ---;t) 
at 

— E Ja 2 a 1 o-a 1 + ^2 K' A2Bl b Bl , (16) 

and 

d 

(ih— -h B2 )b B2 {---O a ---;t) 
at 

= ^J , B 2 A 1 a A 1 +^2,K B2Bl b Bl , (17) 

where h Al — E Al + h v describes the polarization of the 
amorphous network caused by an electron in localized 
state 4> Al mediated by the e-ph coupling, and h B2 = E B . 2 + 
h v . a A2 (• • • a ■ ■ ■ ; t) is the probability amplitude at mo- 
ment t that the electron is in localized state 4> A2 while 
the vibrational state of the atoms is given by normal co- 
ordinates {0 a , ce — 1, 2, • • • , 37V}. &b 2 (- • • & a • • • is the 
probability amplitude at moment t that the electron is in 
extended state £s 2 . 

The e-ph interaction can cause a static displacement 
of atoms. The potential energy shift caused by the dis- 
placement of atoms is: 

av = ^y^kjkXjXk ~y^gpx P , (is) 

jk p 

where g p is the average value of the attractive force 
-dU(r, {Tl n })/dX p ~ Z*e 2 /(4ire e s r 2 ) of electrons act- 
ing on the p th degree of freedom in some electronic state. 
The second term of (TT5)) comes from the e-ph interaction, 
which acts like an external field with strength g p . AV can 
be written as: 

Av = \ E - ~ 4) - \ £ (19) 

jk jk 

where 

^ m^J29p( k ^)mp, m = l,2,--.3/V (20) 

p 

is the static displacement for the m th degree of freedom. 
The constant force g p exerted by the electron on the p th 
vibrational degree of freedom produces a static displace- 
ment for the m th degree of freedom. The deformation 
caused by the static external force of e-ph interaction is 



balanced by the elastic force. A similar result was obtained 
for a continuum model [27] ■ The last term in (|19j) is the po- 
larization energy, a combined contribution from the elastic 
energy and e-ph interaction. 

Owing to the coupling of localized state A\ with the 
vibrations of atoms, the origin of each normal coordinate 
is shifted [35]: 

e a ^e a -ei\ (21) 

where 0* ~ (N A J N)(Z* e 2 / Aire^l J/ {M a u 2 a ) is the 
static displacement in the normal coordinate of the a th 
mode caused by the coupling with localized state A\ , where 
N Al is the number of atoms in region D Al . (j2Tj) leads to 
a modification of the phonon wave function and a change 
in total energy. Using (k~ 1 )jk = (fc -1 )^ and the inverse 
of (fT3|) . one finds that the shift O^ 1 of origin of the a th 
normal coordinate is related to the static displacements 
by: 

0^=^(4-%^°, x° k €D Al . (22) 

k 

The eigenfunctions of h Al are: 

<P N (z) = {2 N N\n 1 ' 2 r 1 ' 2 e- z2 ' 2 H N {z) 1 (23) 

where Hn(z) is the N th Hermite polynomial, 
d a = (MaLUa/hY^Oa is the dimensionless normal coor- 
dinate and Q^ 1 = (M a w CI /?i) 1 ' /2 ^ 1 . The corresponding- 
eigenvalues are: 

4f«> = E° Ai + Y(N* + + £ b Al , (24) 

a 

In an amorphous semiconductor, an electron in state 
A\ polarizes the network and the energy of state |^4i{iV a }) 
is shifted downward by £ b M ~ fc" 1 [Z*e 2 /4ire e s £ 2 ii } 2 . The 
eigenvalues and eigenvectors of h Bl are: 

4f"> =E Bl +Y(N a + l)fru a , (25) 

a 

3JV 

2.3 Static displacement and vibrational amplitude 

In this subsection we compare the relative magnitude of 
static disorder, static displacement of atoms induced by 
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e-ph interaction, and the amplitude of the atomic vibra- 
tions. For a localized state Sa„, one needs to make fol- 



9v 



= d 



l'A 2 



id 



I'JU 



is defined 



lowing substitution in ([T| 
after ©) if p E D A2 and g p = if p £ D A2 . The typi- 
cal value of spring constant k of a bond is k ~ Mu 2 ~ 
Z*e 2 /(4ireo£ s d 3 ), M is the mass of a nucleus, w is a typ- 
ical frequency of the vibrations. The static displacement 
of an atom is x^ ~ g/k ~ (d/^) 2 d. A typical thermal 
vibrational amplitude u 11 is 



u" ~ ^k B T/Muj 2 ~ d(/c B r) 1 / 2 (Z*e 2 /47re e s d) _1/2 . 
The zero point vibrational amplitude is 



^Huj/Muj 2 ~ d(m/M) 1/4 (n 2 /md 2 ) 1/4 (Z*e 2 /47re e s d)^ 1/4 , 

where to is the mass of electron. The e-ph interaction 
for extended states is weak. From both experiments and 
simulations [SDlHTj . the variation of bond length (i.e. static 
disorder) is of order ~ 0.05c?, where d is a typical bond 
length. Now it becomes clear that for amorphous semicon- 
ductors the static disorder is much larger than the static 
displacements of the atoms induced by e-ph interaction. 
The static displacement caused by the e-ph interaction is 
important only when the static displacement is compara- 
ble to or larger than the amplitude of the atomic vibra- 
tions. For weakly polar or non-polar amorphous semicon- 
ductors, the following three statements are satisfied: (1) 
static disorder localizes the carriers in band tails; (2) carri- 
ers in localized tail states have a stronger e-ph interaction 
than the carriers in extended states, and the network is 
polarized by the most localized tail states; and (3) carriers 
in extended states have a weaker e-ph interaction and car- 
riers in extended states are scattered in the processes of 
single-phonon absorption and emission. The small polaron 
theory assumed that e-ph interaction was dominant and 
led to self-trapping of carriers. This assumption is suit- 
able for ionic crystals, molecular crystals and some polar 
amorphous materials. For weakly polar or non-polar amor- 
phous materials, the aforementioned estimations indicate 
that taking carriers to be localized by static disorder is a 
better starting point. 



where Ffj^, -.(t) is the probability amplitude at moment t 
that the electron is in extended state B\ while the vibra- 
tional state of the nuclei is characterized by occupation 
number {N' a , a = 1, 2, • • • , 37V}. 

Substitute Ea.([5S|l and Eq.^ into Eq.® and ap- 
plying J TT t^tf^^ to both sides we obtain 



ih 



dt 



= J2 (MNJIVtMK}) 



■N' 



tt[£ i A ^M )/h 



C#, } (i)e"^ 

(A 2 {N a }\V^ L ph \B 1 {K}) 



B-i—N' 



where 



(A 2 {N a }\Vl r L \A 1 {N' a }) = 



(28) 



(29) 



describes the transition from localized state A\ with phonon 
distribution {• • • N' a • • •} to localized state A 2 with phonon 
distribution {• • • N a • • •} caused by transfer integral Ja 2 a x 
defined in Eq.®. 



{A 2 {N a }\V^ h \B x {N' a }) = 

/n^if a} (E ^v)4^ } 



(30) 



is the transition from an extended state to a localized state 
induced by electron-phonon interaction. 
Similarly from Eq. (|17[) we have 



{N a } 



dt 



= £ (BANaWl^A.iK}) 



2.4 Second quantized representation 

We expand probability amplitude a Al (• • • O a ■ ■ ■ ; t) with 
the eigenfunctions of h Al ■ 



{K} n - lt £ l < } /h 



(26) 



n Ai it{£ 

i N Lr 



+ Y, (B2{N a }\V e E 7 h \B 1 {N' a }) 



F Bi Jt(£ 



{«»} 



where 



(B2{N a }\Vt E \MK}) = 



(31) 



where C^, j (t) is the probability amplitude at moment t 
that the electron is in localized state A\ while the vibra- 
tional state of the nuclei is characterized by occupation 
number {N' a , a = 1, 2, • • • , 37V}. Similarly we expand the 

probability amplitude (• • • a ■ ■ ■ ; t) with eigenfunctions describes the transition from localized state \A% ■ ■ ■ N' a 
of Hbi- to extended state \B 2 ■ ■ ■ N a • • •) caused by transfer inte- 

, gral J' B a , the dependence on {xj} in J' is neglected. 

{B2{N a }\V^ h \B x {N' a }) = 



(32) 

) 
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B 2 B 1 )-^B 1 



(33) 



e) 2 tanh^}x 



is the matrix element of the transition between two ex- 
tended states caused by electron-phonon interaction. Eq. 
and Eq. (|3ip are the evolution equations in second-quantized 
form. The phonon state on the left hand side (LHS) can 
be different from that in the right hand side. In general, 
the occupation number in each mode changes when the 
electron changes its state. 



3 Transition between two localized states 
3.1 J A . as perturbation 



(27T) 



1/2 



£( 



9^) 2 uj 2 a csch 



/3tkOg_ 1-1/2 

2 J 



exp(- 



(35) 



At 'very' high temperature (fc^T > 2.5hu) using tanhx 
x and cscha; « 1/x, (1331) becomes 



W T {A 1 ->• A 



J 2 



In amorphous solids, the transfer integral (0) between two 
localized states is small. Perturbation theory can be used 
to solve Eg. ((28)) to find the probability amplitude. Then where 



XllIzbT 



VLLe 



E^ L /k B T 

A, 



2 



the transition probability |2"3"] from state y to state 



W r (Ai -> A3) = 



^ exp{ -J- [(E° A3 +e+> )- (E A t )] } ) 



AG LL = (E° A3 +£ A >)-(E Ai +£ A i 



b 



(36) 



(37) 



h 2 



1 



/3huj a 

coth }x 



drexp^pft, - {EX +£^)]} 



[exp{^(^ 3 -O^fcsch^cosr^)}-!]. (34) 

We should notice: (i) for localized states we adopt the 
partition ([55)1 for full potential energy (cf. Appendi? jA"]) . 
the LL transition is driven by transfer integral Ja s Ai\ (ii) 
The localized band tail states strongly couple with the 
atomic vibrations 37,38], a carrier in a localized tail state 
introduces static displacements of atoms through e-ph in- 
teraction, so that the occupied localized state couples with 
all vibrational modes. When a carrier moves in or out of 
a localized tail state, the atoms close to this state are 
shifted. In normal coordinate language, this is expressed 
by #q 3 _ ^q 1 m (Ell) for each mode. Thus a LL transi- 
tion is a multi-phonon process; (iii) If we notice that the 
transfer integral J a 3 a x oc e~ R ' 31 ^, where 2/£ = £^ +Ci • 
The product of first two factors in (IMl) is similar to the 
single-phonon transition probability obtained in [7] . In the 
following subsection, we will see that (|3~4"]) is reduced to 
MA theory when reorganization energy is small or two lo- 
calized states are similar: 9 Ai « 9 Al or when temperature 
is high. 



3.2 High temperature limit 

For high temperature {ksT > hTJ, 
to: 



2ttD) 



W t (Ax -»• A 3 ) 



J 



A 3 Aj 

h 2 



-0, 



reduces 



exp^K^+f* 8 )-^ \ LL = 0.2eV, yielding W T ~ 10 12 sec- 



AG LL is the energy difference between two localized 
states. All is the reorganization energy which depends 
on the vibrational configurations {Q A3 } and {0 Al } of the 
two localized states. Because A does not have a deter- 
mined sign for different states and modes, one can only 
roughly estimate 

All ~ k- 1 (Z*e 2 /ATre s e o e) 2 ~ eJ 1 (d/Z) 3 {Z*e 2 /*ireoZ). 

This is consistent with common experience: the longer 
the localization length (the weaker the localization), the 
smaller the reorganization energy. From (|36[) . we know 
that E% l is about 0.01-0.05eV, in agreement with the ob- 
served value[43] for a-Si. ((36)) has the same form as Mar- 
cus type rate ([!]) for electron transfer in a polar solvent 
and in large molecules. Because xq > u v , the vibrational 
energy ku v2 jl is the lower limit of the reorganization en- 
ergy All ~ NAkx^/2. For most LL transitions, All is 
greater than AGll- For less localized states and higher 
temperature, All ~ AGll, then E^ L = A/A + AG L l/2 + 
{AGll) 2 /^A ~ AGll, and the present work reduces to 
MA theory. 

From ab initio simulations[40,4I in various a-Si struc- 
tural models, the distance between the two nearest most 
localized tail states is Ra 3 a ± ~ 3 — 5 A (one or two bond 
lengths). The effective nuclear charge [44] is Z* = 4.29 and 
static dielectric constant[J5] e s = 11.8, Ja z a x ~ 0.02eV 
(Appendix B). The energy dependence AGll between the 
final and initial states affects Wll, mobility and the con- 
tribution to conductivity. For LL transition, the largest 
AGll is the mobility edge D, so that we pick up D/2 = 
0.05eV[12] as a typical AGll- For the most localized tail 
state in a-Si, the localization length [4T)ll4T] is £ » 5A. 
JaxA 2 is estimated in Appendix B. From the force constant [46] 
k ~ dc 4 4, C44 = 81 GPa, a typical reorganization energy 
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If one assumes the same parameters as above, the pre- 
diction of MA theory would be 

(n D /2 or n D/2 + l)J 2 /h{D/2) ~ 10 12 - 10 13 sec ~ x (at 
T=300K), the same order of magnitude as the present 
work, where n D / 2 — {eP D l 2 — is the phonon occupa- 
tion factor. This is why MA theory appears to work for 
higher temperature. 



o 
> 



3.3 Field dependence of conductivity 

For electrons, external field F lowers the barrier of the 
LL transition 5(AGll) = -eF£ - eFR < along the 
direction opposite to the field, where R is the distance 
between centers of two localized states. 

An electric field increases the localization length of a 
localized state. A localized electron is bound by the ex- 
tra force / <~ Z*e 2 u s /4ireoe s d 3 of the disorder potential: 
fu s = h 2 /2m£ i 2 . The relative change <5£ in localization 
length induced by the external field is = —Sf/2f, 

where Sf — —eF is the force exerted by external electric 
field. Thus ~ (eF/2){Z*e 2 u s /[4ire e s d 3 })- 1 > 0. As 
a consequence, reorganization energy A decreases with in- 
creasing F. From A ~ g 2 /k ~ ej 2 Z*e 2 d 3 / (4ire ^), the e -SX/(2X)^ 

one finds: 

relative change S A in reorganization energy A is <5A/A = 
-4<5£/£ w -2eF(Z*e 2 u s /[iireoesd 3 ])- 1 < 0. 

From the expression of E^ L in (|36[) . to first order of 
field, the change in activation energy is 




120 160 

temperature (K) 



200 



Fig. 1. The slope s(T) (see text) for a-Ge: diamonds are 
data[48|. The dashed line is expected from the simple picture 
that the potential energy drops along the inverse direction of 
field for electrons [32 U 50 1. The solid line is calculated from (|39[) . 
To estimate the average over various localized states, we take 
AGll = k B T. 



(1 + a;)- 1 / 2 « l-x/2 w er x l 2 transform (l + ^A/A)" 1 / 2 



s{T) = e( 



Z*e 2 u s , 
4ire e s d 3 (T)' 



R 



Z*e 2 u s 



£(T) v 47re e s ei 3 (T)' 



^d + f). 08) 



ik B T 



(1 



2e{ 



\ 2 {T)' y 4ire Q e s d 3 (T) 



r 



For temperatures lower than the Debye temperature, AGll 
All- It is obvious that SE1 L < 0, activation energy de- 
creases with external field. 

Increasing £ with F also leads to that transfer integral 
J increases with F: since J AsAl oc e~ R31 '^, 
J A3Al (F)JJ A3Al (0) = exp{i? 31 (r 1 -^f 1 )} 
« exp{^ 1 -R3i^/C} > lj where £f is the average local- 
ization length in external field. Using the value of 
J MAl {F)/J MAl {Q) 

= exp{^ 1 i? 3 i(eF/2)(Z*e 2 U V[47reoe s d 3 ])" 1 }- 

For hopping conduction, the conductivity a is esti- 
mated as a = ne 2 /i, for n carrier density and /i = JJ/fc^T 
the mobility, D is the diffusion coefficient of carriers [S]. 
To obtain conductivity, one should average mobility over 
different AGll, density of states and occupation number. 
We approximate this average by AGll ~ ksT. If one only 
considers the contribution from the hopping among near- 
est neighbor localized states, D = R 2 i3Ai Wx(Ai — > A3). 
The force produced by the experimental field is much 
weaker than the extra force produced by the static disor- 
der eF << f, no carrier is delocalized by the external field. 
The carrier density n and the distance R AsAl between 
two localized states are not affected by external field, so 
that a(T,F)/a(T,0) = W T (F)/W T (0). According to (|55jl. 

a(T,F)/a(T,0) = [X(F)/X(0)}-^ 2 [J A;tAl (F)/ J A3Al (0)? 
exp{-l3[E^ L (F)-E^ L {0)}}. Workers often fit experimen- 
tal data in form: a(T, F)/a(T, 0) = exp[s(T)F]. Using 



< 



eTO + R) 
2k B T 



(1 



AG 



LL 



A(T) 



(39) 



According to the percolation theory of the localized- 
delocalized transition [17]) the localization length £ of a 
localized state increases with rising temperature: £(T) = 
£o(l — T/T m )~ 1 (the critical index is between 1/2 and 1; 
we employ 1 here), where T m is the temperature where 
all localized states become delocalized; T m is close to the 
melting point. Then A(T) = g 2 (T)/k{T) = X (l-T/T m ) 4 
and the slope s(T) in exp[s(T)F] increases with decreas- 
ing temperature. Figure [T] is a comparison between the 
observations [15] in a-Ge and the values of present work. 
The parameters used are d = 2.49A, u s /d — 0.1 and 
T m =1210 K. A = 0.2eV is estimated from Z* = 4. Be- 
cause the conductivity comes from various localized states, 
AGll varies from to D. 

The field polarizes the wave functions of occupied states 
and empty states (with a virtual positive charge). A static 
voltage on a sample adds a term to the double- well poten- 
tial between two localized states: 



U(y) 



1 



1 



by - eFy, 



(40) 



where a ~ —k and b ~ fc/xg. To first order in the field, 
the two minima y± and y 2 of (|40[) do not shift. To second 
order in field, the distance between two minima of ([40]) 
decreases by an amount SR — (3b/4a 3 )^/— a/b(eF) 2 . This 
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Fig. 2. Field dependence of conductivity in vanadium oxide. 
The triangle symbols are measured [49] values of a(F)/a(F = 
0) for VO1.83 at 200K. The solid line is the change in carrier 
density n(E)/n(E = 0) computed from Frenkel-Poole model, 
dashed line is the change in mobility fi(F)/fi(F = 0) from 
present work. 



results in a further decrease of reorganization energy [2~T| 
in addition to the direct voltage drop. 

In Figure [2J we compare the observed a(F)/a(F = 0) 
resultsgH] at 200K with the best fit of exp(const-F 1 / 2 ) 
of the Frenkel-Poole model and with exp(s(T)F + dF 2 ) of 
present work (VOi.83 is of special interest for microbolome- 
ter applications [49] ). where 



d = {k B T)- 1 {Z*e 2 /ATre a e s R 2 ){3b/Aa 3 ) v / -a/be 2 . 

We can see from Figure [5] that the change in mobility 
provides a better description of experiments than that of 
the change in carrier density by field in the Frenkel-Poole 
modelEOl. 



3.4 Meyer-Neldel rule 

From the formula for All in the paragraph below ([3"9")) . the 
reorganization energy for LL transitions decreases with 
rising temperature: 

<5A LL (T) = -4ST(1 - T/T m )- 1 T m 1 A LL (T). 

From flU, 

6E^(T) = -25.(0)^(1 -^(l AGIl ^ XLL 



(41) 

If E^ L {T) decreases with rising temperature according to: 



E^(T)=E^(0)(1-T/T MN ), 



(42) 



then the Meyer-Neldel rule is obtained [5T|: comparing (|41l) 
and (|4"2"j) one finds 



MN 



T 

x^ 1 



AG 



LL 



)(1 



AG LL 



^LL 



r 1 . 



According to (|35[) . during a LL transition the vibrational 
configurations of two localized tail states are reorganized. 
A large number of excitations (phonons) is required. A 
temperature-dependent activation energy implies that en- 
tropy must be involved and is an important ingredient 
for activation. The present approach supports the multi- 
excitation entropy theory of Yelon and Movaghar [l2lll3j . 



Table [T] is a comparison of the predicted Meyer-Neldel 
temperature Tmn with observed ones. The number of 
atoms involved in a localized tail state is taken to the sec- 
ond nearest neighbor. The reorganization energy is esti- 
mated from the parameters given in Table [TJ then Tmn is 
estimated from (|4"3"]). Beside NiO, the most localized state 
extends to about l.bd. In NiO the hole of d electron shell is 
localized on one oxygen atom. The localization comes from 
the on-site repulsion in a d-band split by the crystal field. 
The theory agrees well with observations in quite different 
materials. In a typical ionic crystal like ZnO, the localized 
tail states arise from thermal disorder and are confined in 
very small energy range: Ejj ~ ku v2 /2 ~ fc^T ~ 0.025eV 
(T=300K) . The fraction of localized carriers is much less 
than that in an amorphous semiconductor where localiza- 
tion is caused by static disorder. That is why AGll is 
about 10 times smaller than that in an amorphous semi- 
conductor. In ZnO, most of the carriers are better de- 
scribed by large polarons. Taking the melting point as 
the localized-delocalized transition temperature T m is pre- 
sumably an overestimation, so that the computed Tmn is 
too high. 



3.5 Low temperatures 

For low temperature (fc^T < to/10), the argument in the 
last exponential of ([34)) is small. The exponential can be 
expanded in Taylor series and the 'time' integral can be 
completed. Denote: 

f{^) = \{^-e^fcsch^, (44) 



then: 



W T (A 1 ^A 3 ) = r^exp{- 



(43) 



i {3huj a 

exp^E^-^ooth^} 

a 

{J2fMl[S(AG° LL + hiu a ) + 5(AG° LL - fiu a )] 

a 

+ Y J fMf(u a ')^[S(AG° LL + hw a + huj a ,) 

aa' 

+5(AG° LL - hbj a - ftu) a ') 
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Table 1. predicted and observed Meyer-Neldel temperature Ta/jv in several materials (see text). 



material 


T m (K) 




d(A) 


i/d 


AG LL (eV) 


r^theory /-r/-\ 
J-MJV K^J 


T££(K) 


a-Si:H 


1688 


11.9 


2.35 


1.7-1.8 


0.1 


756-890 


499-776 11 


a-Gei_ a; Se2Pb :c 


993 


13 


2.41 


1.7 


0.11 


744 


765 52 


a-(As 2 Se 3 )ioo- ;c (SbSI) :c 


650 


11 


2.4 


1.5 


0.12 


596 


591 53 


ZnO 


2242 


9.9-11 


1.71 


1.5 


0.01 


595 


226-480 54,55 


NiO 


2257 56 


11[57] 


1.75 


1 


0.06 57,58 


1455 


1460-1540 59 



fied. From a linear fit, we deduce T up « fiaJ/2.3fcs, far be- 
yond the more stringent condition T< fEJ/lOks for Taylor 
expansion of the exponential in (|34l) . This is easy to un- 
derstand: low frequency modes are acoustic, the density of 
states decreases quickly with reducing phonon frequency 
(Debye square distribution). T< ftZJ/lOks is derived from 
csch^ < 0.01 for all modes. The exponent of (O is a 
summation over all modes. At T up w fiaJ/2.3fcs, although 
a single phonon seems have higher frequency, because the 
density of states at this frequency is small, the available 

400 800 1, 200 vibration energy is low and VRH is already dominant. 

first peak in phonon spectrum (cm ) 



600- 



e 

o.; 

> 



400- 



200- 



a-Si^Ta-Cu GeSe 



Fig. 3. The present work predicts that the upper temperature 
limit of VRH is proportional to v. The squares are observed 
data in several materials 61 ,62,63,64,65,66,67 , the line is a 
linear fit. 



+ S(AG° LL + hoj a - hoj a ,) + 5{AG° LL - hoj a + hoj Q ,)} H }. 

(45) 

One may say that the transfer integral is reduced by a 
factor exp{ — \J2a(®a 3 ~ ^a 1 ) 2 } due to the strong e-h cou- 
pling of localized states. 

The derivation of (|4"5|) from suggests that below a 
certain temperature T„ p , the thermal vibrations of net- 
work do not have enough energy to adjust the atomic 
static displacements around localized states. The low tem- 
perature LL transition (|45[) becomes the only way to cause 
a transitions between two localized states (it docs not need 
reorganization energy). Since the variable range hopping 
(VRH) 60] is the most probable low temperature LL tran- 
sition, T Mp is the upper limit temperature of VRH. 

On the other hand, the available thermal energy is 
J2j fruijiij, where rij is the occupation number of the j th 
mode. At a low temperature, rij ~ exp(—hujj/kBT), the 
available vibrational energy is determined by the number 
of modes with fko < &bT. The higher the frequency v of 
the first peak in phonon spectrum, the fewer the excited 
phonons. In other words, for a system with higher v, only 
at higher temperature could one have enough vibrational 
energy to enable a transition. Therefore the present work 
predicts that for different materials, the upper tempera- 
ture limit T up of VRH is proportional to the frequency v 
of the first peak in phonon spectrum. 

Figure[3]reports experimental data for a-SijHI], a-Ge[5TJ 
[62] . a-Si02 [63 , 64 , 65 and a-Cu?GeSe3 [6611671 : the upper 
limit temperature T up of VRH vs. the first peak of phonon 
spectrum v. A linear relation between T up and v is satis- 



4 Transition from a localized state to an 
extended state 

The transition probability from localized state A\ to ex- 
tended state Bi is: 

W T {A X B 2 ) = J A^ e -^-EX-sl 1 )/2 
nr 

a 

j\rex P {j(E B2 -E° Ai -£ b Ai )} 
[expi^e^^cos^r} - 1]. (46) 

a 

When a carrier moves out of a localized state, the nearby 
atoms are shifted and the occupation number in all modes 
are changed. Thus a LE transition is a multi-phonon pro- 
cess. 

For 'very' high temperature k B T > 2.5JuJ, (|4"6"[) reduces 

to: 

4 ale 

(47) 

where 

AG° LE = E B2 - (E Ai + £ b A J (48) 

is the energy difference between extended state B 2 and 
localized state A\, and 

n \leKb± 2 ^ 

(49) 
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Here, \le is the reorganization energy for transition from 
localized state A\ to extended state B^. It is interesting 
to notice that the activation energy E% E for LE transition 
can be obtained by assuming As — in All- Transition 
from a localized state to an extended state corresponds 
to that of a particle escaping a barrier along the reaction 
path[36]. 

Formally, the LE transition is similar to LL transi- 
tion. To obtain the former, one makes the substitutions: 

Ja 3Ai -> J' B2Al , et - et e£ and \{e a - + - 

(E A:i + S^ 3 )} -> [E B2 - (E° Ai + £ Al )]. However the physi- 
cal meaning of the two are different. From (|37l) and (|49p. 
we know that Xle is the same order of magnitude as All- 
AG LE is order of mobility edge, which is much larger than 
AG\ L and J'b 2 Ai >> J'a 3 Ai- The spatial displacement of 
the electron in a LE transition is about the linear size 
of the localized state. From (|4"T1) and W£ E becomes 
comparable to W^ L only when temperature is higher than 
k B \E LE - E^HJle/Jll)}- 1 . The mobility edge of 
a-Si is about 0.1 - 0.2eV[40p41]. so that the LL transition 
is dominant in intrinsic a-Si below 580K. However if higher 
localized states close to the mobility edge are occupied due 
to doping, there exist some extended states which satisfy 
AG° LE ~ AG° LL . For these LE transitions, E LE is compa- 
rable to E LL . The LE transition probability is about 10 
times larger than that of LL transition. For these higher 
localized states, using parameters given for LL transition 
in a-Si, W% E ~ lO^sec" 1 . 

According to approximation (i) : Yba < < 1 (Appendix 
A) , (|4"o| is only suitable for localized tail states which are 
far from the mobility edge. ([4T>]) complements Kikuchi's 
idea of 'phonon induced delocalization' [§] : transitions from 
less localized states close to mobility edge to extended 
states. For less localized states, the coupling with atomic 
vibrations is weaker 37, 38], the reorganization energy Alb 
is small. The transition from a less localized state (close 
to mobility edge) to an extended state is thus driven by 
single-phonon emission or absorption [9 , similar to the MA 
theory [7 . Consider a localized state and an extended state, 
both close to the mobility edge. Then AG\ E is small, Wle 
can be large. The inelastic process makes the concept of 
localization meaningless for the states close to the mobil- 
ity edge gHH]. 




'0 100 200 300 400 

Temperature (K) 

Fig. 4. photocurrent time decay index as function of temper- 
ature: circles are data for a-Si:Hp] and a-As 2 Se 3 [l8], dashed 
line is from the transition probability of MA theory, solid 
line is from ()36I47[) in G « A limit. At higher temperature, 
Xle ~ AGle, and the present result reduces to MA theory. 

5 EL transition and EE transition 

If an electron is initially in an extended state |-Bi), the am- 
plitude for a transition can be computed in perturbation 
theory. The probability of the transition from extended 
state \B\) to localized state \A-i) is: 

W T (Bx -> A 2 ) = ^eM~(E A2 +£A 2 -E Bl )} 
cx P {-i £(^ 2 ) 2 coth ^}(h + hi (50) 

a 

where: 

h = J_ dreMj(E A2 + £a 2 - E Bl )} 

1 \ — "V A \2 ftfyjJa 
ex Pl2 Z^(®<x 2 ) csch 2 C0SUJ aT}X 
a 

[E^«bJ 2 + i« Sl ^) 2 coth^l) 

a' 



When a gap-energy pulse is applied to amorphous semi- 
conductors, the transient photocurrent decays with a power- 
law: t r( - T \ The exponent is r(T) = -1 + k B T/E v accord- 
ing to a phenomenological MA type transition probability, 
E v is the Urbach energy of the band tail[T7 | IT5 P ] . @7J) 
leads to [23] r(T) = -3/2 + 2k B T/E v if we follow the rea- 
soning in [1511201132] . Figure H depicts the decay index as 
a function of temperature. At lower temperature the ex- 
perimental data deviates from the prediction of the MA 
theory. At higher temperature and for states close to mo- 
bility edge, A ~ AG and E a ~ AG, the present theory 
reduces to MA theory [201132] . 



CSch COSW a 'T 



a' 

(51) 

and: 

T * = \\Y f ^ K MB 1 Otn 2 ] j^drexp{^(E A2 +£ b A2 -E Bl )}x 
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[exp{^(e) 2 c^^cosc Q r} - 1]. (52) 

a 

For very high temperature fcgT > 2.5fSJ, the EL tran- 
sition probability is 



W EL = u EL e- E ^/ k - T , (53) 



with 



a,' 

*-f L 2M a /w a / 2 8 4 



(^ 2 +gl-^) 2 1 



E 



cscft- 



/3hcu a > 



^« Bl ^ 2 ) 2 coth^ 



-csch 



f3tvjj 



csch — — ]}, (54) 



and 



(55) 



where Xel is order of magnitude of All. Because AG EL < 
0, from the expressions of E EL and -E^ L , we know E^ L 
is smaller than Since K A B Q a f is the same or- 

der magnitude as Ja 3 Ax, the EL transition probability 
is larger than that of LL transition. In a-Si, this yields 
W EL ~ 10 13 — 10 14 sec _1 . When a carrier moves in a lo- 
calized state, the nearby atoms are shifted due to the e-ph 
interaction, so that all modes are affected. Thus an EL 
transition is a multi-phonon process. 

AG° EL < has a deep consequence. From the expres- 
sion ([2~4"[) for £ h A and the order of magnitude of mobility 
edge [12], we know that the energy difference AG EL < 
is order several tenths eV. For extended states with 
\AG el \/\el < 1, we are in the normal regime: the higher 
the energy of an extended state (i.e. AG° EL / Xel more neg- 
ative but still \AG° EL \/ Xel < 1), the smaller the activa- 
tion energy E EL . The higher extended state has a shorter 



lifetime, therefore the time that an electron is able to 
remain in such an extended state is less than the time 
it spends in a lower extended state. An extended state 
with shorter lifetime contributes less to the conductivity. 
For extended states with energies well above the mobil- 
ity edge (such that \AG° el \/Xel > 1), we are in Marcus 
inverted regime (cf. Figure [5]): the higher the energy of 
an extended state, the larger the activation energy. The 
higher extended states have long lifetimes and will con- 
tribute more to conductivity. In the middle of the two 
regimes, AG EL / X E l sa — 1. For these extended states, 
no activation energy is required for the transition to lo- 
calized states. Such extended states will quickly decay to 
the localized states. In experiments, there is indirect evi- 
dence for the existence of this short-lifetime belt. In a crys- 
tal, phonon-assisted non-radiative transitions are slowed 
by the energy- momentum conservation law. In c-Si/SiC>2 
quantum well structure, the photoluminescence lifetime 
is about 1 ms, and is insensitive to the wavelength |70| . 
The photoluminescence lifetime of a-Si/Si02 structure be- 
comes shorter with a decrease in wavelength: 13ns at 550nm 
and 143ns at 750nm[7I]. The trend is consistent with the 
left half of Figure [5] The observed wavelength indicates 
that the energy difference (> 1.66eV) between the hole 
and electron is larger than band gap (1.2eV), so that the 
excited electrons are in extended states. According to ([54")) 
and Figure [5j higher extended states are more quickly 
depleted by the non-radiative transitions than the lower 
ones, so that a photoluminescence signal with higher fre- 
quency has a shorter lifetime. We need to be careful on 
two points: (i) the observed recombination time is order of 
ns, it is the EE transitions that limits EL transition to a 
large extent; (ii) for a quantum well, the number of atoms 
is small, so that the reorganization energy is smaller than 
the bulk. The static displacements may be able to adjust 
at the experiment temperature 2-10K. To really prove the 
existence of the short-lifetime belt, one needs to excite 
electrons into and above the belt with two narrow pulses: 
if the higher energy luminescence lasts longer than the 
lower energy one, the existence of a short-lifetime belt is 
demonstrated. 

In the conduction band, the energy of any localized tail 
state is lower than that of any extended state, so that a 
zero-phonon process is impossible. Because the AG EL < 
0, factor ex-p[\AG EL \/(2kBT)] increases with decreasing 
temperature. On the other hand, other factors in ([rj3"j) de- 
crease with decreasing temperature. Thus there exists an 
optimal temperature T*, at which the transition proba- 
bility is maximum. If one measures the variation of lumi- 
nescence changing with temperature in low temperature 
region, at T* the lifetime of the photoluminescence will be 
shortest. 

A transition from one extended state to another ex- 
tended state is a single-phonon absorption or emission 
process driven by e-ph interaction. The transition prob- 
ability from extended state \B\) to extended state I-B2) 



W{B 1 ^B 2 )= 2 -^^ K b 2 b 1 ) 



M a *uj a > 
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0.1 0.2 0.3 0.4 

E — E (unit: eV) 

E L 

Fig. 5. The non-radiative transition lifetime of extended state 
as function of energy difference Ee — El between initial ex- 
tended state and final localized state. \el = 0.2 eV is esti- 
mated from the data given in Sec. III. When Ee — El = Xel, 
lifetime is a minimum. The short lifetime belt exists in range 
0.1 eV< Ee — El < 0.3eV. The vertical axis is scaled by v^, L . 

~N\ W~, + 1 

l^-5(E B2 -E Bl -ruu a ,) + ^ 5(E B2 -E Bl +tu» a ,)}, 

" (56) 

where N' , = {e^ hu "' — 1) _1 is the average phonon number 
in the a'th mode. In a crystal, (|56[) arises from inelastic 
scattering with phonons. 



5.1 Four transitions and conduction mechanisms 



The characteristics of the four types of transitions are 
summarized in Table [5J The last column gives the order 
of magnitude of the transition probability estimated from 
the parameters of a-Si at T=300K. In a-Si:H, the role of 
hydrogen atoms is to passivate dangling bond, the esti- 
mated rates are roughly applicable to a-Si:H. The rate 
of LL transitions is between two nearest neighbors. For 
an intrinsic or lightly n-doped semiconductor at moderate 
temperature (for a-Si T < 580K, the energy of mobility 
edge), only the lower part of conduction tail is occupied. 
Then AGle is large, and the LE transition probability is 
about two orders of magnitude smaller than that of the 
LL transition. For an intrinsic semiconductor at higher 
temperature or a doped material, AG LE becomes com- 
parable to AG LL , LE transition probability is about ten 
times larger than that of LL transition. The first three 
transitions increase the mobility of an electron, whereas 
EE transition decreases mobility of an electron. 

Although Wel > Wle, the transient decay of pho- 
tocurrent is still observable [18) . The reason is that for ex- 
tended states below the short-lifetime belt, Wel ~ Wee, 
so that an electron in an extended state can be scattered 
into another extended state and continue to contribute to 
conductivity before becoming trapped in some localized 
state. 



If a localized tail state is close to the bottom of the 
conduction band, for another well localized state and an 
extended state close to mobility edge, AGle ~ AGll+D, 
and Wee is one or two orders of magnitude smaller than 
Wel- If a localized state is close to mobility edge, AGle ~ 
AGll, because J' is several times larger than J, Wle 
could be one order of magnitude larger than Wll- The 
probability of EL transition is then one or two orders of 
magnitude larger than that of LL transition. The reason 
is that AG el < 0, E^ L is smaller than E^ L while K'O 
is the same order of magnitude as J. The probability of 
EE transition Wee is about 10 3 times larger than Wll 
(cf. Table [2]) ■ The EE transition deflects the drift motion 
which is along the direction of electric field and reduces 
conductivity. This is in contrast with the LL, LE and EL 
transitions. The relative contribution to conductivity of 
four transitions also depends on the number of carriers in 
extended states and in localized states, which are deter- 
mined by the extent of doping and temperature. At low 
temperature ( k B T < fiZJ/10), the non-diagonal transition 
is still multi-phonon activated process whatever it is LL, 
LE or EL transition. The activation energy is just half of 
the energy difference between final state and initial state 
(cf. ([37]), gll) and ©). 



5.2 Long time and higher order processes 

The perturbation treatments of the four fundamental pro- 
cesses are only suitable for short times, in which the prob- 
ability amplitude of the final state is small. Starting from 
a localized state we only have L — ► E process and L — > L 
process. Starting from an extended state, we only have 
E — >• L process and E — > E process. For long times, higher 
order processes appear. For example L — > E — > L ^> L — > 
E^rL^rE^E^rL etc. Those processes are impor- 
tant in amorphous solids. In a macroscopic sample, there 
are many occupied localized and extended states. If we 
are concerned with the collective behavior of all carriers 
rather than an individual carrier in a long time period, 
the picture of the four transitions works well statistically. 

All four transitions are important to dc conductivity 
and transient photocurrent. In previous phenomenological 
models, the role of LE transition was taken into account 
by parameterizing the MA probability 18,20] . The details 
of LE transition and the polarization of network by the 
localized carriers were ignored. The present work has at- 
tempted to treat the four transitions in a unified way. Our 
approach enhanced previous theories in two aspects: (1) 
the role of polarization is properly taken into account; and 
(2) we found the important role played by the EL transi- 
tion and associated EE transition in dc conductivity and 
in the non-radiative decay of extended states. 



6 Summary 

For amorphous solids, we established the evolution equa- 
tions for localized tail states and extended states in the 
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Table 2. Some features of 4 types of transitions in a-Si. 



transitions 


origin 


phonons needed 


activated 


role in conduction 


probability Wt (sec 1 ) 


L— >L 


J © 


multi 


yes 


direct 


10 iy 


L— >E 


j <nui 


multi 


yes 


direct+indirect 


10 10 - 10 13 


E^L 




multi 


yes 


direct+indirect 


10 13 - 10 14 


E^E 




single 


no 


reduce 


10 13 



presence of lattice vibrations. For short times, perturba- 
tion theory can be used to solve Eqs.© and (0. The 
transition probabilities of LL, LE, EL and EE transitions 
are obtained. The relative rates for different processes and 
the corresponding control parameters are estimated. 

The new results found in this work are summarized 
in the following. At high temperature, any transition in- 
volving well-localized state(s) is a multi-phonon process, 
the transition rate takes form ([TJ. At low temperature, 
variable range hopping appears as the most probable LL 
transitions. 

The field-dependence of the conductivity estimated from 
LL transition is closer to experiments than previous theo- 
ries. The predicted Meyer-Neldel temperature and the lin- 
ear relation between the upper temperature limit of VRH 
and the frequency v of first peak of phonon spectrum are 
consistent with experiments in quite different materials. 

We suggested that there exists a short lifetime belt of 
extended states inside conduction band or valence band. 
These states favor non-radiative transitions by emitting 
several phonons. 

From (|59| one can see that single-phonon LL transi- 
tion appears when states become less localized. In intrin- 
sic or lightly doped amorphous semiconductors, the well 
localized states are the low lying excited states and are 
important for transport. Carriers in these well localized 
tail states polarize the network: any process involving oc- 
cupation changes of well localized tail states must change 
the occupation numbers in many vibrational modes and 
are multi-phonon process. Moving towards to the mobil- 
ity edge, the localization length of a state becomes larger 
and larger. When the static displacements caused by the 
carrier in a less localized state are comparable to the 
vibrational displacements, one can no longer neglected 
the vibrational displacements in The usual electron- 
phonon interaction also plays a role in causing transition 
from a less localized state. In this work, we did not discuss 
this complicated situation. Formally when reorganization 
energy All between two localized states is small and com- 
parable to the typical energy difference AGll, the present 
multi-phonon LL transition probability reduces to single- 
photon MA theory. 

The multi-phonon LE transition discussed in this work 
is for well localized states, it is a supplement to the the- 
ory of phonon-induced derealization which is concerned 
with the less localized states close to the mobility edge. 
As discussed in Appendix A, when the atomic static dis- 
placements caused by a carrier in a less localized state is 
comparable to the vibrational amplitude, LE transition 



could be caused either by a phonon in resonance with the 
initial and final states or by the transfer integral Jb 2 a x ■ 

We thank the Army Research Office for support under MURI 
W91NF-06-2-0026, and the National Science Foundation for 
support under grants DMR 0903225. 



A Approximations used to derive evolution 
equations 

Usually in the zero order approximation of crystals (spe- 
cially in the theory of metals), the full potential energy 
En U a (r - R„ - u£) is replaced by V c = J2 a U n (r - R n ), 
where R n and u™ are the equilibrium position and vi- 
brational displacemet of the nth atom. One then diago- 
nalizes h c — — fi, 2 V 2 /2m + V c , such that all eigenstates 
arc orthogonal. Electron-phonon interaction Y^=i x j Jx~ 
slightly modifies the eigenstates and eigenvalues of h c or 
causes scattering between eigenstates. This procedure works 
if e-ph interaction does not fundamentally change the na- 
ture of Bloch states. As long as static disorder is not strong 
enough to cause localization, Bloch states are still proper 
zero order states. However one should be careful when 
dealing with the effect of static disorder {u™}. It must 
be treated as perturbation along with e-ph interaction. If 
we put static disorder in potential energy and diagonalize 



-ft 2 V 2 /2m + £ n t/ n 



Ft,, 



u"), the scattering 



effect of static disorder disappear in the disorder-dressed 
eigenstates. The physical properties caused by static dis- 
order e.g. resistivity is not easy to display in an intu- 
itive kinetic consideration based on Boltzmann-like equa- 
tion: eigenstates of h s are not affected by static disorder 
{u™}. By contrast, computing transport coefficients with 
eigen states of h s is not a problem in Kubo formula or its 
improvment{34 . 

We face a dilemma in amorphous semiconductors. On 
one hand the static disorder {u™} is so strong that some 
band tail states are localized, static disorder must be taken 
into account at zero order i.e. diagonalize h s ; on the other 
hand from kinetic point-of-view the carriers in extended 
states are scattered by the static disorder which should 
be displayed explicitly E n u " ' dU n (r — R n )/<9R n rather 
than included in the exact eigenstates of h s . 

The very different strengths of the e-ph interaction in 
localized states and in extended states also requires differ- 
ent partitions of the potential energy ^ n U n (r — lZ n — u"), 
where lZ n = R n + u" is the static position of the n th 
atom in an amorphous solid. Molecular dynamics (MD) 
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simulations [37 ,38] show that the eigenvalues of localized 
states are strongly modified (about several tenth eV) by 
e-ph interaction, while the eigenvalues of extended states 
do not fluctuate much. It seems reasonable that for lo- 
calized states we should include e-ph interaction at zero 
order, and put it in the zero-order single-particle poten- 
tial energy (just like small polaron theory [221I23] ) . while 
for extended states e-ph interaction acts like a perturba- 
tion (similar to the inelastic scattering of electrons caused 
by e-ph interaction in metals). 

The traditional partition of full potential energy is 

]T^-^ n ,u n ) = ]>>(r-ft n ) + 5> — . (57) 

n n— 1 j—1 3 

In this ansatz, static disorder is included at zero order. 
Localized states and extended states are eigen states of 
h s . The second term of ([5"7j) . the e-ph interaction, is the 
unique residual perturbation to eigenstates of h s . It causes 
transitions among the eigen states of h s i.e. LL, LE, EL 
and EE transitions, to lowest order transition is driven 
by single-phonon absorption or emission. Since all attrac- 
tions due to static atoms are included in h s , two types of 
transfer integrals Ja 2 a x (from a localized state to another 
localized state) and J'b 2 a x (fr° m a localized state to an 
extended state) do not exist. One can still use Kubo for- 
mula or subsequent development 34 to compute transport 
coefficients. 

However partition ([5T|) obscures the construction of a 
Boltzmann-like picture for electronic conduction where 
various agitation and obstacle mechanisms are explicitly 
exposed. The elastic scattering caused by static disorder 
is hidden in the eigen states of h s . To explicitly illustrate 
the elastic scattering produced by static disorder, one has 
to further resolve the first term of ([5"T]) into 

M jV 
£t/(r-ft n ) = ^t/(r-R n ) 

n=l n=l 

+ J2 u? ■ dU n (r - R„)/c9R n . (58) 

n 

([57]) is also inconvenient for localized states. The e-ph in- 
teraction for a carrier in a localized state is much stronger 
than in an extended state [571155] . It is reflected in two as- 
pects: (i) a localized carrier polarizes network and pro- 
duces static displacements for the atoms in which the 
localized state spread; (ii) the wave functions and cor- 
responding eigenvalues of localized states are obviously 
changed (the change in eigenvalues can be clearly seen in 
MD trajectory [371138]). To describe these two effects, in 
perturbation theory one has to calculate e-ph interaction 
to infinite order. 

Taking different partitions for localized states and ex- 
tended states is a practical ansatz.Yov a localized state, we 
separate 

Y i U{T-K u ,x%)= Y, U(r-TZ n ,u^) 

n ngDjij 



+ Y tf(r-ftp,u?)> (59) 

P<£D Al 

where Da ± is the distorted region where localized tail state 
(j)A 1 spreads. For a nucleus outside Da ± , its effect on lo- 
calized state A\ dies away with the distance between the 
nucleus and Da % - The second term leads to two trans- 
fer integrals Ja 2 a ± (induces LL transition) and J' B A (in- 
duces LE transition) in the evolution equations of localized 
states. 

Since for a well-localized state 4>a 1 , the wave func- 
tion is only spread on the atoms in a limited spatial re- 
gion Da ± , {<f>Ai \ ' \<f>Ai ) = 0, where • stands for the sec- 
ond term in the RHS of ([ST)]) . In calculating Ja 2 a x = 
(<j>A 2 \ ■ \<t>A t ) and J'b 2Ai = (Cflal • |<AAi), it is legitimate 
to neglect u p , ■ » J2 p <£d a U{r — Tip)- The change in 
potential energy induced by the atomic vibrational dis- 
placements is fully included in the first term in RHS of 
([5T)1) . Because the wave function cf>A 1 of localized state A\ 
is confined in Da x , one can view cf>A 1 as the eigenfunc- 
tion of h° Ai = ~h 2 V 2 /2m + £ n6l , U(r - K a , <) with 
eigenvalue Ea 1 ({u™, n 6 Z?^ }). The rest of atoms outside 
-Dai act as boundary of 4>Ai_- A carrier in localized state 
4>Ax propagates in region Da 1 and is reflected back at the 
boundary of Da 1 .Since static disorder is fully contained in 
h Ai , in the present ansatz, localized carriers are free from 
clastic scattering of static disorder. 

For a less localized state §a x , its wave function spreads 
over a wider spatial region Da 1 ■ With shift to the mobil- 
ity edge, J2neD Al u " ' VJ7(r — 1Z n ) become smaller and 
smaller, eventually comparable to X) p ^_d a U(r — 7t p ) and 
J2 p <£d a u v ' ^U(r — 7?. p ,u p ). For carriers on these less 
localized states, their polarization of the network is weak, 
and the atomic static displacements are comparable to the 
vibrational amplitudes. Entering or leaving a less localized 
state does not require configuration reorganization, and 
the reorganization energy becomes same order of mag- 
nitude as vibrational energy. For such a situation, only 
a phonon in resonance with the initial and final states 
contributes to the transition. The multi-phonon processes 
gradually become the single-phonon processes, although 
the driving force is still the transfer integral induced by 
Y^ p etD A — 7?.p)- O ne obtains the MA theory. 

For a less localized state, if we treat J2neD A u " ' 
VC/(r-K n ) ) Ep^ 1 ^-^ P )andE p ^ 1 <-VC/(r- 
7?p,uP) in the same foot as perturbation, the phonon- 
induced derealization naturally appears and is accom- 
pany with EL transitions induced by transfer integral Jb 2 a x 

For extended states, to construct a kinetic description, 
([571) is suitable partition of full potential energy. The sec- 
ond term of ([571) . the e-ph interaction, causes EL and EE 
transitions in the evolution equation of extended state. 
The elastic scattering of the carriers in extended states 
induced by static disorder can be taken into account by 
two methods: (i) Using the eigenvalues and eigenfunctions 
of h s in Kubo formula; (2) if we wish to deal static disorder 
more explicitly, we can apply coherent potential approxi- 
mation to ([55]) . In this work we do not discuss these issues 
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and only concentrate on the transitions involving atomic 
vibrations. 

In deriving (|6I7[) . we neglected the dependence of ex- 
tended state on the vibrational displacements, then 
Vj^b-l = 0. Since v n /v e ~ 10~ 3 (where v e and v n are typ- 
ical velocities of the electron and nucleus), m/M ~ 10~ 4 
(m and M are mass of electron and of a typical nucleus) 
and x/d ~ 1CP 2 — 1CP 1 (x and c? are typical vibrational 
displacement of atom and bond length), one can show 

E^MXv^xv^j or EjC^MjOv^ « 

E P ^ Al I/(r-«p,uP)0 Al or EjXiPU/dXAfa. There- 

fore terms including Vj^Ai or V^^i can be neglected. 

To further simplify the evolution equations, we need 
two connected technical assumptions: (i) the overlap inte- 
gral Yb 2 Ax between extended state £b 2 and localized tail 
state 4>a 1 satisfies Yb 2 a x ~ N^/N « 1 and (ii) over- 
lap integral Sa 2 a x between two localized tail states sat- 
isfies Sa 2 a x << 1- Assumption (i) means that we do not 
consider the localized tail states very close to the mobil- 
ity edge and consider only the most localized tail states. 
Condition (ii) is satisfied for two localized states which do 
not overlap. It means we exclude the indirect contribution 
to conductivity from the transitions between two localized 
tail states with overlapping spatial regions. 

For two localized states A 2 and A\, Sa 2 A! « 1 if Da 1 
and Da 2 do not overlap. The terms multiplied by Sa 2 a 1 
can be neglected for localized tail states which their spatial 
regions do not overlap. What is more, the transfer integral 
is important only when the atoms p ^ Da ± fall into Da 2 
or Ai = A 2 , 



A J __ A r-. 



Y a Ai J A 2 A l +Wa 2 CIA 2 



(60) 



where Wa 2 = J d 3 r\^A 2 1 2 E p {d A2 U{r - Tip, uP) only af- 
fects the self energy of a localized state through a,A 2 ■ Com- 
paring with Ea 1 and with h v , Wa 2 can be neglected [55]. 



B Scales of the coupling parameters of four 
transitions 

For these most localized tail states, the wave functions 
take the form of c/)a 1 ~ e~l r_TC - 4 i Z^ 1 . Ja 2 a x is estimated 
to be -(N A2 Z*e 2 /4ire e s £)(l + R 12 /0e~ Rl2/ ^ where av- 
erage localization length £ is defined by 2£ _1 = ^-f 1 
i?i2 is the average distance between two localized states, 
e s is the static dielectric constant, Z* is the effective nu- 
clear charge of atom, Na 2 is the number of atoms inside 
region Da 2 ■ Later we neglect the dependence of Ja 2 a x on 
the vibrational displacements {x} of the atoms and con- 
sider Ja 2 Ax as a function of the distance R± 2 between two 
localized states, localization length £i of state 4>a x and 
localization length £2 of state <j>A 2 ■ If a localized state is 
not very close to the mobility edge, its localization length 



£ is small. The overlap between it and an extended state 
Yb 2 Ax ~ Nai/N ma y be neglected. 

If we approximate extended states as plane waves £bj ~ 
e ik »i r , then K' MBi ~ (Z* e 2 U /47re e^ 2 )(l - ike^)' 1 . 
So that K' A2Bi u/ 'J 'a 2 a x ~ e Rl2 ^u/t;, where u ~ 

y/ksT '/Mui 2 or yj h/Mcj is typical amplitude of vibra- 
tion at high or low temperature. The distance between 
two nearest localized states is ~ several A in a-Si, and 
K'a 2 Bi u i s several times smaller than Ja 2 Ax - If we again 



-ikt 



2 . J' B2Ai is Of 



approximate extended state as plane wave £ B 

J'b 2Ai ~ (Z*e 2 /4ne s s U 1 )(l + ik B2 U 1 ) 
the same order of magnitude as Ja 2 a ± - J'b 2 a 1 does not 
create transitions from an extended state to a localized 
state. The asymmetries in (fTQ)) and (jlip come from the 
different separations (I5T)1) and (1551) of the single particle 
potential energy for localized states and extended states. 
One should not confuse this with the usual symmetry 
between transition probabilities for forward process and 
backward process computed by the first order perturba- 
tion theory, where two processes are coupled by the same 
interaction. If we approximate extended states and 
£b 2 by plane waves with wave vector k\ and k 2 , Kb 2 Bi ~ 
Z*e 2 UK 3 /(47reo£ s [K+«(fc 2 -fci)]), where k ~ (e 2 /e )(dn/dfi) 
is the Thomas-Fermi screening wave vector. In a lightly 
doped or intrinsic semiconductor, n is hundreds or even 
thousands times smaller than 1/a, a is bond length. Since 
for most localized state, localization length £ is several 
times a, therefore Kb 2 b ± ~ ( k £) 2 ^_a 2 -Bi > ^ s mucn weaker 
than three other coupling constants. 



C LE transition at low temperature 

The LE transition probability for low temperature ksT < 
fuU/ 10 can be worked out as in (|43|) . Denote: 



(61) 



and the result is 



W T (A 1 -> B 2 ) 



2n,J' 



12 

B 2 A 



1 exp{- 



-f3AG 



LE ~ 



1 /3huJ a 

exp{--]T(e) 2 coth— }x 

a. 

^hE{w a )\[&<A&LE + + S(AG° LE - fiu a )] 

a. 

+ Y fLE{u a )fLE(Ua')^[${AG° LE + ftuJ a + tiU) a >) 
aa' 

+5(AG LE — huj a — hu> a i) 



-S(AG° LE + hjj a - hcu a ,) + 6(AG LE - huj a + thu a /)] H }. 

(62) 
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D EL transition at low temperature 

For low temperature (k B T < fiu/10), one can expand the 
exponentials in (l5Tj) and (|52"j) into power series. Then the 
integrals can be carried out term by term. To 2-phonon 
processes, the transition probability from extended state 
\Bi) to localized state \A 2 ) is 

W T {Bi -> A 2 ) = ^eM-^(E°A 2 +£ h A 2 E Bl )} 



-5{E° A2 + S\ - E Bl + huj a „ - hu a ,„) 



-5{E° A2 + £ b M - E Bl - hw a n + nuaf,.)- 



6(E° A2 + £ b - E Bl - hu a „ - Hu a n,)] 



{ K A 2 B X f CJch 2l 3huJa ' 
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{ 2^ [ 2M a >u a > 4 



• coth ]csch- 

2 J 2 



[6{E° A + £ b M - E Bl +Hu a ,) + S(E° A2 + E\ 2 - E Bl - hoj a , )] 



where 



[6(E A2 +£ b A2 -E Bl +2fiw a ,) 
5(E° A2 +£ b A2 -E Bl -2huj a ,)} + ---}, (63) 



f E M = \{6tfcsch^. (64) 



E 



{Kt Bi et?f 



(1 - csch 2 —^ )^2f E L(u a ») 



[8{E\ 2 + £\ 2 - E Bl +hu al ,) + 8{E\ 2 + E\ 2 - E Bl - hu a 



la qA 2 \2 



1 \^^ K A2B 1 ) h , (kzbo: 



4 ^ 1 2M Q «w Q - 

a'a" 



■ coth 



CSCh JELy-^a") 



X [ S ( E °A 2 + £a 2 - E B 1 + ?U*J a ' + tiU) a » 



+S(E° A2 + £\ - E Bl + hw a , - fiua») 



+S(E° A + £ b - E Bl - fUslaf + feda") 



+S(e° A2 + e\ - e Bi - hu a , - hue.)] 



(.4 



{2-csch z ^Y L -)]x 



E fEL(u a ")f E L{0J a "') 



x [s(e° A2 + e\ - e Bi + hoj a „ + hw 
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